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Extra dimensions are a common feature of beyond the Standard Model physics. 
In a braneworld scenario, local physics on the brane can depend strongly on the 
brane's location within the bulk. Generically, the relevant properties of the bulk 
manifold for the physics on/of the brane are neither local nor global, but depend on 
the structure of finite regions of the bulk, even for locally homogeneous and isotropic 
bulk geometries. In a recent work p], various mechanisms (in a braneworld context) 
were considered to stabilize the location of a brane within bulk spaces of non-trivial 
topology. In this work we elaborate on and generalize that work by considering 
additional bulk and brane dimensionalities as well as different boundary conditions 
on the bulk scalar field that provides a Casimir force on the brane, providing further 
insight on this effect. 

In D = 2 + 1 (D = 5 + 1) we consider both local and global contributions to the 
effective potential of a 1-brane (4-brane) wrapped around both the 2-dimensional 
hyperbolic horn and Euclidean cone, which are used as toy models of an extra- 
dimensional manifold. We calculate the total energy due to brane tension and elastic 
energy (extrinsic curvature) as well as that due to the Casimir energy of a bulk scalar 
satisfying both Dirchlet and Neumann boundary conditions on the brane. In some 
cases stable minima of the potential are found that result from the competition of at 
least two of the contributions. Generically, any one of these effects may be sufficient 
when the bulk space has less symmetry than the manifolds considered here. We 
highlight the importance of the Casimir effect for the purpose of brane stabilization. 
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I. INTRODUCTION 

Several outstanding problems remain in fundamental physics. To name a few: unifica- 
tion of the gauge interactions of the Standard Model with the gravitational interaction, the 
enormous hierarchy between the observed gravitational and weak scales, M pl /M weak w 10 17 , 
as well as the hierarchy among fundamental fermion masses (e.g. m T /m e w 3500), and 
the smallness of an apparent cosmological constant, A, where y/X ~ H m 10 _33 eV~ 
10- 44 M wcak W. 

Extra spatial dimensions might exist and may help to explain some, if not all of these 
issues: they are an essential feature of string theory; "large" extra dimensions have been 
employed in attempts to explain the weak (e.g. [2H1]) and flavor (e.g. [5]) hierarchies; there 



[1] As usual we have set H = c = 1. 
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have also been attempts to address the dark energy problem via infra-red modifications of 
gravity that employ extra dimensions (e.g. [SHE])- 

Common experience and laboratory experiments tell us we don't appear to live in more 
than three spatial dimensions; if more exist they should be hidden in some way. One 
explanation might be that the Standard Model fields are confined to submanifolds of lower 
dimension or perhaps the extra dimensions are compactified (i.e. the manifold has non- 
trivial topology) to length scales that have remained inaccessible by experiments. The actual 
situation may be a combination of both. Furthermore, since topology is a global property 
of the space-time manifold it cannot be determined by the local Einstein's equations - 
interestingly, though, topology can have a large influence on local physical processes within 
the manifold. 

In what follows below we will consider scenarios where the full space-time manifold is the 
direct product of 4D Minkowksi space-time, Ai A , (or similarly, FRW) with a d-dimensional 
extra manifold, Q d , but certainly the discussion could extend to other scenarios. With d — 1 
only the topology needs to be specified, as the geometry is trivial. From experience we know 
geometry can play an important role as well, therefore we begin our discussion with d = 2. 

Let us first consider a cylinder with topology (R x S 1 ). It can be obtained from Euclidean 
2-space (E 2 ) by modding out by a group (r) of 1-dimensional discrete translations. The 
local geometry remains the same as the covering space (i.e. homogeneous and isotropic), 
however the presence of fields indicates that some (global) symmetry has been broken. For 
example, the components of field momenta are quantized in the compact dimension, while 
they remain continuous in the infinite dimension. That is, local experiments are sensitive 
to the global structure of the manifold. Even so, no special places exist on this space, as the 
results of such experiments are insensitive to where on the manifold they are performed. 

This should be contrasted with a the two-dimemsional hyperbolic horn ("H 2 /r), obtained 
from hyperbolic 2-space H 2 by modding out by the same T (see Figure [I]). Fields on this 
space indicate that both rotational and translational symmetry are broken; field modes with 
non-zero excitation in the compact dimension are highly inhomogeneous in z and are doubly 
exponentially suppressed beyond the point on the manifold where the "circumference" of 
the compact dimension becomes smaller than the mode wavelength. In other words, modes 
of a certain momentum can only be sent a finite distance down the horn. Since translation 
invariance is broken, there exists a notion of absolute position and the results of local 



FIG. 1: Horn (top) and Cone (bottom) topologies. 



experiments will vary along the length of the horn. The main lesson here is that even if the 
local geometry of a manifold is homogeneous and isotropic, the physics on it will generally 
not be. 

The physically relevant scale in this example is the horn circumference, a quantity that 
is sensitive to non-local characteristics of the manifold, yet is insensitive to the whether or 
not the manifold is actually entirely a horn that extends to infinity with infinite volume. 
For this reason we term such a quantity regional, a quantity that could not have been 
calculated simply using local geometric information, yet did not require complete knowledge 
of the manifold's structure. The hyperbolic horn is of great utility because it has a lot of 
the salient features of more generic manifolds, in particular the "cuspy" regions of compact 
hyperbolic manifolds (CHMs)' 2 l It therefore seems that what may appear as a global effect 
can actually be attributed to the regional structure of the manifold one considers. 

Another simple but non-trivial topology is that of a cone, a two-dimensional Euclidean 
space with a "wedge" removed from it. It is obtained by identifying 9 <H- 9 + 2tt(1 — 5), 
when working in polar coordinates, r and 9 (see Figure [T]). Topologically, the space is given 
by M 2 — {0}, i.e. the point at the tip of the cone is removed from the manifold. Again, 
translation invariance is broken, as the distance from the tip of the cone is imprinted on the 

[2] CHMs provide an appealing geometric solution for the hierarchy problem (see [9] ) ; they could be considered 
a hyperbolic, d > 1, version of models discussed in [2] or [4], wherein all but one modulus is fixed. CHM's 
have also found use in string theory (see e.g. [TP]). 
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modes. The cone offers a complementary example to the horn as it has zero curvature and 
is a better approximation to regions of manifolds that are geometrically flat and end in a 
vertex. 

We would like to focus our attention on braneworld scenarios and the local physical 
consequences of regional properties of the bulk manifold. They would affect, e.g. vacuum 
energy of a brane itself, as it is generically proportional to the inverse of the brane volume, 
which can vary throughout the manifold. Furthermore, the amplitude for brane fields to 
interact with fields in the bulk (e.g. the gravitational field) would also be sensitive to the 
position of the brane, as the structure of bulk modes are generally inhomogeneous. Finally, 
as we shall explore in this work, there are local forces on branes arising from e.g. the Casimir 
effect of bulk fields interacting with the brane. 

Before exploring any interesting phenomenological prospects, however, we must ask what 
determines the brane's location within the bulk and what ensures it is stabilized. From an 
effective 4D perspective the brane's position within the bulk appears as a series of massless 
scalar fields, one for each of the brane's codimenensions. Observational constraints on fifth 
forces tell us that such fields should have a mass > 10~ 3 eV [H] if they couple to matter with 
gravitational strength. In any case, stability is a reasonable requirement of any system and 
thus the brane's effective potential should have a stable minimum in which it sits. Stabilizing 
a brane in this context should be contrasted with other mechanisms (e.g. [12]) that have 
been used in the context of stabilizing the moduli of the bulk - here we are concerned with 
determining and stabilizing the brane position within a bulk that is assumed to be stable. 

In this work we focus on different local (geometric) and global contributions to the brane 
potential and find that they depend on the manifold's non-trivial structure in different ways. 
From classical intuition, we expect that a brane's non-zero tension will provide an energy 
proportional to the brane volume. Furthermore, there will generally be an elastic energy 
associated with the brane's extrinsic curvature and how it couples to the intrinsic bulk 
curvature. 

Due to local interactions with the brane, inhomogeneities in bulk field modes will be 
induced - this will be approximated here by an effective boundary condition on the bulk 
field(s) at the location of the brane. Therefore if global translation invariance of the bulk 
manifold is broken, the vacuum energy of bulk fields (a global quantity) will depend on 
the location of the brane, providing a location-dependent force, i.e. the Casimir effect [13]. 
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Likewise, the Casimir energy of Standard Model fields living on the brane can play a role 
if the brane's topology is non-trivial. Calculation of Casimir energies is subtle because it 
strongly depends on the bulk and boundary geometry, topology, dimensionality, field type, 
and boundary conditions. The magnitude of the effect can be estimated on dimensional 
grounds; however in order to get the overall sign and precise magnitude of the energy, a full 
calculation usually must be performed. Having explicit analytic expressions for the fields 
modes (as one has on the horn and cone) makes this task much more tractable. 

Though using Casimir energies to stabilize bulk moduli is not a new idea (see e.g [HHIB]), 
to our knowledge they have not been used to stabilize the location of a single brane within 
the bulk nor have the regional properties of manifolds been exploited for their use. We would 
especially like to emphasize its role because, in the case of a 3-brane embedded in a higher- 
dimensional bulk, the brane geometry is trivial, as are the corresponding energies. The 
Casimir effect may then be the only mechanism available to stabilize the brane's location in 
this case. 

In summary, the gradients of the total energy result in a net force, which can stabilize 
the location of a brane. In this work, we generalize the results of pQ, confining our attention 
to a codimension-1 brane wrapped around either the 2-dimensional hyperbolic horn or the 
Euclidean cone, but extending our analysis to additionally include the Neumann boundary 
condition, as well as analogous calculations in D = 2 + 1 (i.e without the Minkowski dimen- 
sions). On both geometries (in D = 5 + 1) we suppose that all Standard Model fields are 
confined to a codimension-1 brane as pictured in Figure [2] or [7| but that they are free to 
propagate in the compact (i.e. "universal" [32]) dimension. This is possible if the "circum- 
ference" of the brane is small enough so that the first excited "Kaluza Klein" mode has not 
yet been made accessible by experiment. 

We have calculated the total energy of these systems as the sum of tension and elastic 
energies (i.e. those local to the brane and associated with its geometry) as well as the 
(global) vacuum energy from a bulk scalar field, 0, that satisfies a boundary condition on 
the brane. We also parametrize the vacuum energy of Standard Model fields confined on 
the brane. For each geometry, we first consider total space-time dimension D = 2 + 1 then 
consider the addition of three Minkowski dimensions in D = 5 + 1. While the 2+1 D cases 
are interesting for academic reasons, they also illuminate how the different contributions 
to the brane potential scale differently with the number of spatial dimensions as well as 
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provide examples of the how the explicit poles of the vacuum energy depend significantly 
on the even/oddness of D. 

In section[TT]we analyze the hyperbolic horn, considering both dimensionalities and bound- 
ary conditions. There we discover, in both dimensionalities, that a competition between 
tension and bulk Casimir energies can result in a stable brane position, however only for 



the Dirichlet case. In Section III we perform the analogous calculations for the cone. A 
stable minimum is possible for all boundary conditions in D = 2 + 1 and can result from the 
competition between the tension and elastic energies (i.e. purely from geometric effects), 
but is also ensured by the Casimir energy from the scalar. In D = 5 + 1 a global minimum 
occurs for the Dirichlet case while only a local minimum is possible for Neumann. In Section 



IV we conclude with a summary of results and discussion. 



II. HYPERBOLIC HORN 



A. Preliminaries 



1. The Model, Energy Contributions, and Bulk <p Solutions 

Here we will consider a general spacetime dimension D — m + 2 + 1, where m indicates 
the number of Minkowski spatial dimensions. That is, the the full spacetime manifold is 
j\Z[m+i x ^2 jY _ ^he [j ne element can be written in coordinates such that 

ds 2 = r l { ™ +1) dx»dx u + e- 2z/z *zld6 2 + dz 2 , 

where we identify (6) <H- (8 + 2ir), and we have chosen our coordiantes such that the 
horn circumference at z = is 2irz*. Furthermore, from here on we shall work in units 
where the intrinsic length scale of the horn, z* = 1. The model is illustrated in Figure 
[2| suppressing any Minkowksi spatial dimensions, where the codimension-1 brane wraps 
around the compact direction and resides at coordinate Zb- In D = 5 + 1 (m = 3) we assume 
the Standard Model fields to be confined to the brane, free to propagate in the compact 
direction. We will now enumerate different contributions to the brane's effective potential 
on the horn. 
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FIG. 2: A partial embedding diagram of the horn. The codimension-1 brane is pictured with 
coordinate z b . The manifold is infinite in extent in the ±z directions. 



Energy Contributions 

There will be a contribution from a non-zero brane tension, a, that provides an energy 

E tcn = / d 1+m x^/W\a = 2ne- Zb V M a , (1) 

where 7^ is the induced spatial metric on the brane, and the volume of the regulated 
Minkowski spatial slice is denoted by Vm- There will also be an "elastic" energy contribution 



due to the extrinsic curvature of the brane, K ab (see Appendix E2a). We write this as an 
expansion in scalars built from K ab , namely 



E CUIV = j d 1+m x^\ (hK 2 + h 2 K ab K ab + ...)= 27ie- Zb V M (h x + h 2 + . 



(2) 



where the hi are parameters describing the energy cost for deforming the brane within the 
manifold, and . . . represents other possible scalars^ 3 ' that are higher order in K ab . Note that 
the effect of a and the hi are completely indistinguishable, at least as far as the total energy 
of the system is concerned. Without loss of generality, then, we shall set hi — and encode 
all geometric effects in a. It is already clear that these effects behave monotonically with z b , 
and are thus insufficient to stabilize the brane by themselves, but as we will see, the Casimir 
energies from the brane or bulk can provide a complementary opposing force. 

The Casimir (or vacuum) energy in the bulk scalar on a constant-t hypersurface is 

E = I V^dE{0\T^\0) (3) 



[3] In principle we should also include possible coupling between the extrinsic- and Riemann curvature tensors, 
however in this geometry those terms appear with the same Zf,-dependence, thus at this level they are 
indistinguishable from the terms with only extrinsic curvature. 
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from which one obtains the canonical result Eq = | Xli^i' where i is a general mode index. 
It is well known that this sum is infinite; to extract the physically relevant energy we employ 
the zeta-function regularization technique (see e.g. [TS]) and write 

i 

where /i is a renormalization scale. The generalized sum is finite for large enough s, which 
is analytically continued to zero. In some systems the terms in the energy that determine 
the Casimir force remain finite even in the limit s — > 0, while in others they diverge and 
must be dealt with explicitly; in either case it is ultimately the renormalization of geometric 
parameters in the full theory (including gravity) that remedies the situation. 

In the D = 5 + 1 model there will also be a vacuum energy associated with the Standard 
Model fields living on the brane, sensitive to the size of the brane in the compact dimension. 
Though we will not calculate this effect explicitly, by dimensional considerations we expect 
its associated 4D energy density to scale approximately inversely proportional to the fourth 
power of this length scale, namely 

KSM / rN 

Po,sm - 4 (5) 

(27re Zb ) 

where ksm is a dimensionless coefficient (recall that we are working in units where z* = 1). 



Bulk (f) Solutions 

The action for a real massless scalar field in D space-time dimensions is 

S 4> = \J d D x^V,<pVy (6) 
whose variation with respect to <fi yields the Klein-Gordon equation on this space-time, 

□0 = (-81 + V 2 ) = (7) 

where 
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and Vx is the Minkowski-space Laplacian. The positive frequency modes, Mi, are 

m = Ae-^-^-^Z^z) , (9) 
where n£Z and A\ is a normalization constant. Z n ^ k satisfies 

KkW - Z' n ^z) + (u, 2 -p 2 - n 2 e 2 ^) Z n>k (z) = , (10) 

where p = |p| is the momentum in the Minkowski directions. We find that the boundary 
conditions and normalizability imply a real k > 0, defined by the dispersion relations 

\/p 2 + k 2 if n = 

- = { 7 (ii) 

' p 2 + k 2 + \ otherwise 



To make the problem more tractable, we regulate the infinite spatial volume of the horn 
by truncating the space at z = Zl, where zl < £&, and impose there a Dirichlet boundary 
condition on the neld' 4 l At the appropriate place in the calculation we take the limit 
Zl — > —oo to recover the full horn spacetime. For the n ^ modes we find (see Appendix 



Al) 



e z/2 [I-ik (\n\e ZL ) I ik (\n\e z ) - I ik (\n\e ZL ) I_ ik (\n\e z )} when z L < z < z b . 
e z l 2 K ik (|n|e 2 ) when z b < z . 



(12) 

The n = modes do not contribute to the Casimir force, as we show explicitly in Appendix 



B3; however this can be understood from the following argument. From (10) we see that 
these solutions behave as e z l 2 sin k(z — zl,r), and the energy density per unit z in these 
modes is proportional to ^/\g\ZQ k = sin 2 k(z — Zl,r), where zl and zr are simply there to 
temporarily regulate the size of the space in the ^-direction. This is the same form as the 
energy density of scalar modes in a 1-dimensional box; as the box size goes to infinity we 
know there is no Casimir force due to translation invariance. Thus the n = modes on 
the horn do not contribute to the Casimir force, essentially, because they don't "feel" that 
translation invariance of the manifold is broken. 



[4] This makes the modes discrete and normalizable. The choice of Dirchlet boundary condition is arbitrary 
and should have no bearing on the end result. 
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For all n, the spectrum of k are determined by the boundary conditions at z&: 

Z n , k (z b ) (Dirichlet) , 
or 

Z' n ^(z b ) (Neumann) . 
Finally, the field may be expanded in terms of the solutions as 



(13) 



E 



t * 

a\U{ + a-M; 



(14) 



where a\ and a\ are the creation and annihilation operators of modes labeled by the set of 
quantum numbers, i. As usual, the vacuum state is defined as a\ |0) =0 and the operators 



satisfy the commutation relations, [a\, Oj] = = [a], at], and [ai, aj] = <5ij. The positive- 
energy eigenfunctions are normalized using the Klein-Gordon norm (assuming a discrete 
spectrum) : 



,ti 



(15) 



where £ is a spacelike hypersurface, and n M is a unit timelike vector normal to it. 



2. Procedure for Calculation of Casimir Energy on the Horn 



Contour Integral Representation of Sum 

As we generally have no explicit expression for the {k}, we will calculate Q in part using 
a contour integral representation, as elucidated in [18] . To this end, the following will be of 
use for both the D = 2+1 and 5+1 models (i.e. m = and m = 3): 

E (* 2 + i) ^" = hl dk (? 2 + i) I ln A »« < 16 > 

{k} 

where 7 is a counter-clockwise contour that encloses the entire spectrum (the positive real 
axis) and A n (k) are appropriate mode-generating functions whose roots correspond to our 
spectrum of k (for details see Appendix [Bj or e.g. references [TS] or [IS]). With properly 
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chosen generating functions, one can show that the contour may be deformed to give 

E (* + i) ^ - f~ *(*-^± - a..) a, 

{fc} 

We have derived a set of mode generating functions for n ^ and n = separately 
(Appendix [B]), and have shown that the n = modes are independent of and so we omit 
n = /rom an?/ sums we encounter. We find the mode generating functions for both sides 
of Zb to effectively combine to give the log of the total generating functions 

{In [I k (\n\e Zb ) K k (Inle**)] (Dirichlet), 
(18) 
In [I' k (\n\e Zb ) K' k (\n\e z »)] (Neumann), 

where any ^-independent terms have been omitted as they do not contribute to the Casimir 
force. 



General Analysis of Eo(s) 



The divergent piece from Eq(s) needs to be separated out so that its physically relevant 
part may be revealed through an analytic continuation in s, else the divergent quantities 
need to be explicitly absorbed through some renormalization of parameters in the full theory. 
Because the summand is even in n and n = isn't counted, we may take 

oo 

E- 2 E a*) 

n n=l 

so that in general we have 

„ (m + l) 

r°° ( \\ ^i~ s f) 
E (s) oc J2 J i/2 dk [k 2 - 4 J ^ ln ^n{ik) (20) 



Given (18), an analytic continuation of (20) in s is impossible in its current form, so we 



approach it using a uniform asymptotic expansion, closely following the procedure in [18J. 
We know the divergences occur at large k and n, but the idea is to isolate their asymptotic 
behavior, taking them to infinity simultaneously while keeping their ratio fixed. If this 
asymptotic behavior can be understood analytically, it is then straightforward to handle the 
divergent contributions so that the remaining finite part may be calculated numerically. To 
this end, we decompose Eq(s) into a sum of its asymptotic and finite part: 

E (s) = E™(s) + E*? (21) 
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We change variables, defining 



and 



Our sum/integral therefore becomes 



y = 



(22) 



(23) 



E (s) cx^xr 1 " 



n=l 



(2x 6 )-l 



Axi 



(ro+1) 



9_ 



In A n (iyx b ) . 



(24) 



From the uniform asymptotic behavior of the two modified Bessel functions (see e.g. [2H]). 
we define the asymptotic part of the log of the generating function as 



ln[I yXb (x b ) K yXb (x b )} = In 



2x by /y 2 + 1 



2x h 



m+2 j, , 

+ Y (Dirichlet) 

pi (y^y 

m+2 i _n 

9j{V ) /AT \ 

- > -V r— Neumann 



(25) 
(26) 



where the fj(z) and <?j(z) are rational functions of z, given in Appendix [C] We find it most 
convenient to capture the divergences of (24) by defining' 5 ^ 



poo / -I 



(m+l) , 



0J/ 



lnf/^ (x b ) K yXb (x b )}. (27) 



for the Dirichlet case, and similarly for Neumann. We will see that the integral may be 
performed analytically, however the sum must be done numerically and so an analytic con- 
tinuation again becomes challenging. We will then find it convenient to further decompose 
Eq s (s) into divergent and "remainder" parts, Eq"(s) and E™" 1 . In summary, the total energy 
(or energy density) of the system will be given by 



E = E tCQ + Et(s) + E™ + E^ (D = 2 + 1) (28) 

P = Ptcn + Po iv (*) + Po° m + PT + Po, sm (D = 5 + 1) (29) 



[5] This turns out to be easier than doing an expansion in inverse powers of x b of the entire integrand. 
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B. D = 2 + l 



Here there are no Minkowski dimensions, so p = 0. From Q and (17) we have 



En(s) 



2.s 



- COS TCS 



°° />oo / 1 

E / . * ( e - A 



1/2-8 



d_ 
dk 



In Ajfc(zfe) 



(30) 



where In A n (ik) is defined by (18) in both the Dirichlet and Neumann cases. 



1. Dirichlet Condition 



Many of the details of this section are included in Appendix [D] and maybe be easily 



generalized to the other cases. Considering (25) and (27), we define the asymptotic part of 
the energy as 



2s 00 poo / 1 \ 1/2-s o _ 

E as ( S ) = -^cosvr^x^ / dy(y 2 -— 2 ) -\n[I yXb (x b ) K yXb (x b )} (31) 



7T 



„=1 - r (2* 6 )- 1 \ 

After analytically continuing towards s = we find 



E as ( S ) 



(2/i) 2s ^ 1 + 16x 2 b + (49 - 10s)x i b + 64^ 

COS ITS 2^ 



4 



1 + Ax 



2^/2+s 



(32) 



n=l V x ~ ^b) 

As oc n, the divergences here will come from the terms in the summand that go asymp- 
totically as Xb or a;^ 1 ; we extract these by defining 



(33) 



where we have expanded in inverse powers of x b to obtain 

9 -26s 



-fi 2s COS TTS 



n=l v 



128 b 



e Zb 9 , /13 

1 e~ Zb 

24 128 



1 



9" " 7 " 2s ~ ln 



and (evaluated with s = 0) 



7-ircm 
E 



E 

n=l 



1 + 16x 2 + 49x£ + 64xg 
4(l+4x 2 ) 5/2 



^ 1 2 X6 + U8 Xb 



(34) 



(35) 
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Asymptotically, we find 

{_7+||ln2 x e -z b ag z _^ _ 

256 ^ 
"llb^ 3 ) X e ~ 3Zb aS Z b ->■ +°° • 

where we have used a lowest order Euler-Maclauren expansion in the first line, and a Taylor 
expansion in x^ 1 in the second. After setting s = 0, integrating by parts' 6 ', and returning 
to the variables k and Xb one arrives at 

E o n = ~ E / dk ~, 771/2" ( ln [ J * ( x b) R k fa)] - Hh (x b ) K k (x b )} ) (37) 



Asymptotically, we find 



2.4 x 10 3 x e Zb as Zb — > — oo , 
^o n ~ . (38) 

— 1.4 x 10 3 x e 326 as z b — > +oo . 

We have learned that £"o lv contains an explicit 1/s pole which is proportional to e~ Zb ; we 
have also verified this using a complementary heat kernel analysis (see Appendix [E]). The 
presence of the pole and the arbitrary scale \i is indicative of the fact that the Casimir energy 
cannot be measured in isolation; it is only a contribution to a total energy that is itself finite 
[21] . One must identify a parameter that is to be renormalized - since this offending term is 
proportional to the brane length, a renormalization of the brane tension will suffice, with \x 
setting the renormalization scale. In so doing, we choose to go beyond minimal subtraction 
and absorb finite corrections into the tension as well, hence we find the asymptotic behaviors 



(see Appendix D) 



2irz+e Zfc / Z * x a ren as Zh — > — oo , 
J5(*)~< (39) 

where we have restored the horn length scale, z*. We find that the total energy is unbounded 
from above as z& — > +oo and likewise for ^ -oo as long as cr ren > 0, and therefore a 
global minimum exists. E(zb) is plotted in Figure [3] 



[6] This is advisable for numerical purposes. 
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FIG. 3: In D = 2 + 1 with a Dirichlet boundary condition on (j), various potential contributions as 
a function of Zh- For illustration we have chosen a ren = 0.02 [z+ 2 ]. 



2. Neumann Condition 



Here we repeat the analysis of the previous section, only now using the Neumann speci- 



fication in (18). Altogether, we find the asymptotic behaviors 



2irz*e Zb / Z * x <7 ren as — > — oo , 
E(z b ) ~ { (40) 

We find that the Casimir energy is unbounded from below for zj, — > +oo (the sign of the 
energy is flipped with respect to the Dirichlet case, in this limit), and brane tension cannot 
compensate for this, therefore no stabilizing mechanism exists in this case. E(zb) is plotted 
in Figure [4] without brane tension. 



C. D = 5 + l 



Here the previous analysis is repeated, adding three Minkowksi spatial dimensions, i.e. 
m = 3. For ease of normalization we compactify these dimensions using a torus of fun- 

1/3 

damental side length, Vy[ , which is taken to infinity to recover a true global Minkowski 
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FIG. 4: In D = 2 + 1 with a Neumann boundary condition, the total energy as function of zj,. Here 
we have chosen cr ren = 0. 

space. In this limit, the Minkowski momentum, p, becomes continuous. The sum over the 
components of p turns into an integral as 

V = ^ Ap^ Vm_ f = Vm f 2 { ) 

From Q we now have the effective 4D vacuum energy 

2s 00 /• / i \ 1/2— s 

^(^) = ^EE/^ 2 K + fc2 + i ( 42 ) 

n=l {k} J V ' 

We compute the p-integral by first performing a Mellin transform, 
which allows us to directly integrate over p, yielding 

, , a 2 * r(»-2) /„ i\ 2 -' 
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We now represent the sum over {k} as a contour integral using (17), giving 



Po [s) 



pL s — 2)sin7rs 



/i 2s 



2-s 



d_ 
dk 



In A n (ik) 



32tt 2 



1+s G" 7 -* ( " /2) ))S/ /2 *r"i 



2-s 



d_ 
dk 



lnA n (ik) (45) 



where ip is the Digamma function, and we have omitted terms higher order in s as they 
vanish in the s — > limit. 



1. Dirichlet Condition 



Considering (25) and (27) we define 



,,2s / /q \ \ 00 /-oo / i \2-« 

" / 1 i „ / ° „, - 1 /„\ \ \ ™4-2s / t , / „,2 1 



„2 



d ~ 

x gyWyx, (x b ) K yXb {x b )\ (46) 

The integral may be performed analytically, yielding a cumbersome expression which we omit 
for brevity. Here the divergences come from terms in the summand that behave asymptot- 
ically as xf,xf or x®. We isolate those terms in (having already performed the sum on 
n) 

Pa v (s) = — (-— + 7" " -21n7rue- z " + ^(-V2)l - ^C'(~2)e 2 ^ + ^C'(-4)e 426 
Ho KJ 2048tt 2 V 12 s P ^ V '') 64vr 2SV ; 32tt 2S v ; 

(47) 

From this we numerically obtain 

= j™ - P^OO) (48) 

which asymptotically behaves as 

{—1.5 x 10~ 4 x e~ Zb as 2 b — > — oo , 
(49) 
6.8 x 10 5 x e 2zb as z b -)■ +oo . 

Lastly, we return to the variables k, x b and integrate by parts to give 
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FIG. 5: For D = 5 + 1 with the Dirichlet condition on (f), the total 4D energy density in units of the 
horn curvature scale (z+), as a function of the brane position, z b - For illustration, we have taken 
the brane tension, cr re n = 

0.09 [z~% and KSM = 0. 
Asymptotically, one may show 

{-2.0 x 10 -5 x e~* b as z b -> -oo , 
(51) 
4.4 x 10 5 x e 2zb as z b — > +oo . 

In contrast to the 2+1 case, the pole in Eq 1v (s) is a constant and is simply discarded. 
Altogether we have the asymptotic behaviors 

p(z b ) ~ { (52) 
(|7r 2 C'(-4) + K SM ) x e Azb ' z * as z b -)• +oo . 

where finite quantum corrections have been absorbed by <7 ren , and l/27r 2 ('(— 4) 0.04. As 
in the D = 2 + 1 case, here there is a global minimum, assuming a ren > and assuming 
the value of k sm doesn't alter these results (see Figure [H]). For a = 0(1) [z~ 5 ], the effective 
mass for the brane's position modulus is 0(z~ l ). 
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FIG. 6: For D = 5 + 1 with the Neumann condition on cj), the total 4D energy density in units of 
the horn curvature scale (<£*), as a function of the brane position, Zf,. Here we have chosen <7 ren = 0, 
and ksm = 0. 



2. Neumann Condition 



Here we proceed as in the previous section, only using (26). Altogether we find the 
asymptotic behaviors 

{2irz ir e~ Zb ^ z * x cr ren as z b — > — oo , 

(53) 
leh* (-^ 2 C'(-4) + k sm ) x e 4 ^* as z b -> +oo . 

Once again, the sign of the vacuum energy in <p is nipped with respect to the Dirichlet 
case in the Zf, — > oo limit. It is thus unbounded from below, and no value of brane tension 
can create a stable minimum, however a sufficiently large and positive Ksm would accomplish 
this (see Figure [6]). 
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FIG. 7: A partial embedding diagram of the cone. The brane is pictured at coordinate r b . As 
opposed to the horn, this manifold is geometrically flat and is infinite in extent only in one direction. 

III. EUCLIDEAN CONE 
A. Preliminaries 

1. The Model, Energy Contributions, and Bulk <j) Solutions 

Once more we work in general spacetime dimension D — m + 2 + 1, where m indicates 
the number of Minkowski spatial dimensions. That is, the topology of the full spacetime 
manifold is Ai m+1 x (IR 2 — {0}). The line element can be written in coordinates such that 

ds 2 = 7 1 { p +l) dx^dx v + dr 2 + r 2 d9 2 

where we identify 9 <H- 9 + 2n (1 — 5), i.e. the conical space is defined by the deficit angle 
2tt8. Again, a 1-brane (4-brane) is wrapped around this space in D = 2 + 1 (.0 = 5 + 1) 
and resides at coordinate distance r& away form the cone tip, as illustrated in Figure [7j As 
opposed to the horn, the cone is flat, having no intrinsic length scale associated with it. 

Energy Contributions 

Here, a non-zero brane tension (a) contributes an energy given by 

E ten = J d 1+m x^/W\o- = 2tt (1 - 5) r b V M o- (54) 
where 7^ is the induced spatial metric on the brane, while the elastic energy characterized 



by the extrinsic curvature of the brane, K a b (see Appendix E2b) is 



O curv = / d 1+m x^/W\ (hK 2 + h 2 K ab K ab + ...) = 2?r(1 % M (ht+h 2 + .. .) (55) 

J 'b 
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where the hi are parameters describing the energy cost for deforming the brane within the 
manifold, and . . . represents other possible scalars'^ that are built from terms higher order 
in K a b. Note that the effect of each of these terms is completely indistinguishable from 
one another, at least as far as the total energy of the system is concerned. Without loss 
of generality, then, we shall define h% + hi + ■ ■ - = h. As opposed to the horn, the tension 
and elastic energies appear with markedly different dependences on the brane location. 
Interestingly, it is already apparent that these two energies are sufficient to stabilize the 
brane under the physically reasonable assumption that a, h > 0. 

The zero-point energy in the bulk is again zeta-function regularized and is given by Q. 
As we will see below, a check on our results is available in D = 2 + 1 when 5 = as this is 
simply the Casimir effect of the 1-sphere. In D = 5 + 1 we parametrize the vacuum energy 
of brane- dwelling Standard Model fields, from dimensional considerations, by 

KsM (Ka\ 

Po, sm ^ — 74 (56) 

(2n(l-6)r b y 

where ksm is a dimensionless coefficient. 



Bulk Solutions 



Again, the massless bulk scalar field obeys the Klein-Gordon equation, 

□0 = (-d 2 t + V 2 ) = (57) 

where 

B 2 1 8 1 <9- 
V 2 = V 2 + — + -— + —— (58) 
x dr 2 r dr r 2 d6 2 



and V 2 is the Minkowski-space Laplacian. The positive frequency modes, u- 1} are found to 



be 



Ui = A i e-< ut -»*'- vff )R lhk (r) 

[7] Coupling to the Riemann curvature tensor is trivial as it is zero for this geometry. 



(59) 
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where v = n/(l — 5), n G Z, Ai is a normalization constant, and i? ni fc satisfies Bessel's 
equation 

K,k(r) + \K k {r) + (k 2 - £\ R n>k (r) = (60) 

whose solutions may be written as a linear combination of the Bessel functions of the first 
and second kind, J u (kr) and N u (kr), or the Hankel functions H^\kr) and H^\kr). The 
positive frequency dispersion relation is 



vV + k 2 (61) 



where p = |p| is the momentum in the Minkowski directions. Note that while u doesn't 
have an explicit n-dependence, the quantized k's do depend on n non-trivially. 

To make the problem more tractable, we regulate the infinite spatial volume of the cone 
by truncating the space at an outer coordinate r out ^> r b , and impose there a Dirichlet 
boundary condition on the field. At the appropriate place in the calculation we take the 
limit r out — > +oo to recover the full spacetime. 

For the interior solutions we use Bessel functions and find that for n ^ the only ones 
that are normalizable are 

Rn^k{r) = Ju{kr) (0<r<r b ) (62) 



while, for n = 0we are allowed (See Appendix A 2) 



Rn=o,k{r) = c x J v (kr) + c 2 N u (kr) (0 < r < r b ), (63) 

From here on we set c 2 = 0, however we find no physical argument as to why this must 
necessarily be true' 8 l For the exterior solutions we use Hankel functions for later convenience, 
and find (for all n) 

R n , k (r) = (kr out ) (At) - F« (Ar out ) (kr) (r b < r < r out ), (64) 



[8] In fact, we must only ensure that the Hamiltonian is Hermitian and that the modes are normalizable; 
there is nothing inherently wrong with singular modes unless, of course, physical quantities diverge. The 
ambiguity in the value of C2 represents an ignorance of what lies at the cone's tip; presumably it would 
be determined uniquely if the structural details of the tip were resolved. 
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where we have already imposed the Dirichlet condition at r out . The spectrum of k are 
determined by the boundary conditions at r b : 



= < 



R n ,k(r b ) (Dirichlet) , 
R' nk (rb) (Neumann). 



(65) 



2. Procedure for Calculation of Casimir Energy on the Cone 



Contour Integral Representation of Sums 

Again, we have no explicit expression for the {k}, and so we calculate (|4]) in part using 
the contour integral representation 

yV- 2 *+™ = _ cos7r(f -s) r°° dk k i-2s +m d lnAn{ik) (66) 
{k} vr J dk 

where the contour has already been deformed. Lastly, we must specify the mode generating 



functions, which are derived analogously to the case of the horn (see Appendix B4). We 
find these functions for both sides of r b to effectively combine to give the log of the total 
generating function 

\\n[I u (kr b )K u (kr b )} (Dirichlet), 
lnA n (ik)=l (67) 

I In [I I (kr b ) K' v (kr b )] (Neumann) . 

where any ^-independent terms have been omitted as they play no role in the Casimir 
effect. 



General Analysis of E (s) 

We need to separate out the divergent parts from E (s) so that they may be either 
analytically continued in s or explicitly absorbed through some renormalization(s). Because 
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the summand is even in n, we will have contributions to the energy of the form 

■lnA (iA:) (68) 



^ n=0) oc 



oo o 

dk k l ~ 2s+m 



and 



E, 



(n?40) 



OC 



E 

n=l 



8k 



oo o 
dfc k l ~ 2S+m 



dk 



In A n {ik) 



(69) 



Given (67), an analytic continuation of these expressions is impossible in its current form. 



We want to decompose E (s) into a sum of its asymptotic and finite part: 



E (s)=E™(s) + E$« 



(70) 



For n = this simply means performing an expansion of In Ao(ik) in inverse powers of 
k. Following [22J, we use the large k expansion to define 

m+l 



In [K (kr b )I (kr b )} = \n 



In [K' {kr b )I' {kr b )}=\n 



2kr b 



2kr b 



+ E^ 

m+l D 

+ E Bj 



(Dirichlet) 



(Neumann) 



(71) 



(72) 



{kr b y 

where the first few coefficents are A 2 = 1/8, A* = 13/64, B 2 = — 3/8, B 4 = —27/64 and 
Aj = Bj = for all odd j. We thus capture all the divergences by taking 



g* (n=0) ^ 



oo o 

ok 



In [if (Ar 6 )/ (A:r 6 )] 



(73) 



for the Dirichlet case, and similarly for Neumann. For n ^ 0, we employ a uniform asymp 
totic expansion (as done for the horn) by first making the variable change 

k 



y 



V 



therefore our sum/integral becomes 

oo 



■2s+m 



n=l 



8 

dy y'-^-^lnA^iuy) 



(74) 



(75) 



From the uniform asymptotic behavior of the two modified Bessel functions (see e.g. [20J) 
we define the asymptotic part of the log of the generating function to be 

m+2 j, , \ 

fjiyn) 



In [K v (uyr b ) I v (uyr b )\ = In 
In [K' v (uyr b ) I' u (uyr b )} = In 



1 yVrg + 1 
2v y 2 r 2 



+E- 

m+2 

+ E 

J'=l 



7/J 



(Dirichlet) (76) 
(Neumann), (77) 
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where the fj(z) and gj(z) are given in Appendix [Cj We therefore capture all the divergences 



of (75) by taking 

00 roc q ^ 

K S ( ^ 0) oc V^- 2s+m / dy y^+^ln [K„ (vyr b ) I v {vyr b )\ . (78) 

, 1 Jo °y 



n=l 

pfin 



for the Dirichlet case, and similarly for Neumann. In all cases, Eq is obtained from the 
difference of E (s) and Eq s (s) To recap, we have in total 

E = E tcn + E curv + E* s (s) + £ fin (D = 2 + 1) (79) 

P = Ptcn + Pcurv + Po(s) + + P 0, SM (£> = 5 + l) (80) 

B. D = 2 + l 



Here there are no Minkowski dimensions, so p = and m — 0. From Q and (66) we 
find 

2s / poo o oo „ \ 

E (s) = -— cosvrs / dA; k 1 ' 23 — kiA {ik) + 2 V / dk k 1 - 23 — In A n (ifc) (81) 
27r afc n=1 Jo ok J 

1. Dirichlet Condition 



The n = modes: In this case we isolate the divergent part as 

^a S (n=o) = _P_ CQS7rs dk k^—ln [I (kr b ) K (kr b )] (82) 
27r Ja ir ok 

which vanishes after taking both Air — > and s — > 0, the order being irrelevant. The finite 
piece is then found numerically (after integration by parts) 



E?> Ji f°°dk 

^ Jo 



In [/ (kr b ) K (kr b )} - In [I (Ax b ) K (Ax 6 ))] 



1 4 x 10" 2 

(83) 

r 6 



The n^O modes: Here we have (with = z/y) 

.2s °° /-oc a ^ 

^ Jo 9j/ 



1 ' ' " ' ^3 + 7 + ^ + In \fir b (1 - 5)] ) (85) 



24(1- 5) r b 128r fe V 2s 
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The finite piece is 

1 r°° r -i 

El a (r ^ 0) = +-Y j dk In [I v (At 6 ) tf„ (At 6 )] - \n[I v (kr b ) K v (kr b )} (86) 
= ^fi (87) 

n 

where Af(S) is a numerical coefficient; we find jV(0) ~ 9.7 x 10~ 4 (when the cone is opened 
fully) and that its magnitude decreases to zero'^ monotonically as 5 — > 1 (in the limit the 
cone closes in on itself). 

Summary: A 1/s pole remains whose existence we verify using a heat kernel analy- 
sis in Appendix [E] and, based on how this term scales, we identify h as the parameter 
that must be renormalized. Taking the total energy as the sum of the tension energy, the 
renormalized curvature, and the vacuum contributions we find 

m , on n , 2vr(l-*)ft rcn 1.3XHT 2 1 1 (l-ijln/in 

E(n) « ar (i - <) no + — + Yijr ^ 7b (88) 

where we have also absorbed finite quantum corrections in h ren . The scale \x remains in the 
argument of the logarithm, as it must for it be unitless. This scale, however, is arbitrary and 
completely degenerate with the measured value of h ren , and may be chosen to be anything 
convenient. Then, without loss of generality, we define [i as the scale at which h ren = 0. 
The logarithmic term dominates in both the r b — > and r b — > oo limits and, because of its 
overall minus sign, the energy has a global minimum (see Figure [8]). Having a > simply 
steepens the potential, however it isn't required for stability. 

Finally, as a check on our results, we set a = 5 = 0, to recover the scalar Casimir energy 
of a 1-sphere with a Dirichlet boundary condition (e.g. [18, 23J) 

fi^Eoifx- 1 ) w 6.8 x 10" 4 (89) 



[9] This is consistent with the fact that in this limit ln[-] becomes a better approximation to ln[-] since it is 
defined as an expansion in large v. 
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FIG. 8: For D = 2 + 1 with the Dirichlet boundary condition on cf>, the total energy in units (fi) 
where h ren = 0, as a function of the boundary position, r b . For illustration, we have taken 5 = 1/4 
and the brane tension, a = 10~ 12 [/x 2 ]. 



2. Neumann Condition 



Here we proceed as in the previous section, using (67). We find the total energy to be 



r,/ n n « 2tt (1 - S) h Ten 0.25 1 1 5(1-5), 

E(r b ) *2ir(l-5)r b a+ —± ^77"^ vox — 90 

n r b 24 (1 - 5) r b 128r 6 

As in the Dirichlet case, the log term dominates in both limits of r&, meaning there is 
again a global minimum' 10 ' in the potential (see Figure [9j). As a check on our results, we set 
a — 8 = 0, to recover the scalar Casimir energy of a 1-sphere with a Neumann boundary 
condition [El [23] 

^Eodi' 1 ) w -0.18 (91) 



[10] One may notice a large difference in scales between the Dirichlet and Neumann case, namely that the 
minima lie at around 10 4 /j, _1 and 10 _4 /i _1 , respectively. This is traced back to the logarithmic behavior 
of the potential, therefore the minimum depends exponentially on 0(1) numbers that are either > or 
< 0. 
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FIG. 9: For D = 2 + 1 with the Neumann boundary condition, the total energy in units (//) where 
hren = 0, as a function of the boundary position, r^. For illustration, we have taken 5 = 1/4 and 
the brane tension, a = 10 5 [/x 2 ]. 

C. D = 5 + l 

Here, we treat the addition of three Minkowksi spatial dimensions (m = 3) in the same 



manner as Section |II C[ regulating the size of the Minkowksi spatial slice by compactifying 

M 



1 /3 

these dimensions using a torus of fundamental side length, VJ . Using the same logic and 



the result (66) we find the effective 4D vacuum energy 

^TTn r(g ~ 2 \ sin7r g y f°° dk k A ~ 2s ^- In A n {ik) 
16tt 5 / 2 r(s - 1/2) ^ J dk nK ' 



Po[s 

1071 

,2* 



64tt 2 V \2 



ir " ( \ ■ A t --<■■( L / 2 )) ) >< 



i o °<J poo o 

dk fc 4 " 2s — \nA (ik) + 2 V / rffc fc 4 " 2s — lnA n (ifc) 



(92) 



31 



1. Dirichlet Condition 



The n = modes: Here we define the analytic part as 

,2s / /g 



as (n=0) 

Po 



9 



rfA; A; 4 ~ 2s --In [J (Ax 6 ) ^ (*r 6 )] 



13 



2048vrV 4 V s 



1 15 
+ 2 In fir b + — - 7 - ■?/>( ~ V 2 ) 
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(93) 



where in the second line we have chosen' 11 ^ Atr = r 6 x . The finite piece is then 
1 



fin (n=0) 
Po 



Air a 
4 ^ 



. . , , } dk A; — In [J (kr b ) K (kr b )] 



+ 



L 



4 ^ 



dA; k 



dk 



In [J (Ax 6 ) K (kr b )\ - In [J (kr b ) K (kr b )} 



-5.4 x 1(T 4 x 



(94) 



The n^O modes: Here, with k = uy we have 

2s / /g 



as (n^O) 

Po 



n=l ^ 



13 



1 



2048.V ,7---21n[2,(l-^] + ^(-./ 2 ) + 



(9 ~ 
199 



7C'(-2) 



C'(-4) 



64vr 2 (1 - <5) 2 r 4 32tt 2 (1 - 5f r\ 
Leaving us with a finite part given by 

j oo 



fin (n^O) 

Po 



8tt 2 



?1=1 



dA; A; d 



In [J, (Ax 6 ) tf„ (Ax 6 )] - ln[/„ (kr b ) K v (kr b )} 



AT (6) x 



(95) 



(96) 
(97) 



where Af(5) is a numerical coefficient. It is about —2.8 x 10 4 for 5 = and its magnitude 
decreases to zero monotonically as 5 increases towards 1. 



[11] In the D = 2 + 1 cases we were able to simply analytically continue Air — > 0. Here we must choose it 
to be finite (but arbitrary) and Air = r^ 1 is convenient for displaying our results. To be clear, the sum 

_ ^as (n-0) _|_ ^fin (n-0) . g j nc l e p enc l en t Q f j n t nc g _j. Q nIm t. 



32 




FIG. 10: For the D = 5 + 1 Dirichlet condition, the total 4D energy density vs. brane position 
using an arbitrary unit of length, L, as a function of the boundary position, r*,. For illustration, 
we have taken S = 1/4, a = 0.05 [L~% h = 0.01 [L~ 3 ], k sm = 0. 



Summary: Taking the total energy as the sum of the tension, renormalized curvature, and 
vacuum energy contributions, we have 



, s n „ rN 2tt (1 - 8) h Tea 229 13 In [2tt(1 - 8)} 

p{r b ) « 2tt 1 - <J r 6 <7 + — ^ ^ + 2 4 1 \ 4 h 98 

r;, 8192^^ WzAir^r^ 



7C'(-2) , 1 fl v „ 



+ C'(-4) + ksm 



64vr 2 (l -5) 2 r^ 16vr 4 (1 - 5) 4 r\ \1 ' 

where l/27r 2 C'(— 4) » 0.04. Here the 1/s poles have canceled when the n = and n 7^ 
modes were added together. As r& — >• the energy is unbounded from above for all 
5, neglecting k$m- Therefore, if a positive brane tension exists the energy has a stable 
minimum. Since the cone has no intrinsic length scale we choose use an arbitrary scale, L, 



as a reference with which to plot our results in Fig 10 Interestingly, even in the absence of 
the Casimir effects of the bulk scalar, a stable minimum is achieved if we take both h and 
a to be positive, a physically reasonable assumption. 
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FIG. 11: For the D = 5 + 1 Neumann condition, the total 4D energy density using an arbitrary unit 
of length, L, as a function of the boundary position, r b . For illustration, we have chosen 5 = 1/4, 
a = 0.05 [L-% h = 0.1 [L~ 3 ], k sm = 0. 

2. Neumann Condition 



Summary: He we proceed as in the previous section, only using (67). The total 4D energy 
density is found to be 



p(r b ) ps 2n(l -S)r b a + 
13C(-2) 



2n(l-6)h 



7309 



n 



1228807r 2 r^ 



+ 



27 In [2tt (1 - 5)} 
10247r 2 r^ 



(99) 



+ 



+ 



irX '(-4) + ksm 



64tt 2 (1 -5) 2 r^ 16vr 4 (1 - 5) 4 V 2' 

Here the energy is unbounded from below as r b — > for all 5, neglecting ksa/. A global 
minimum apparently doesn't exist when ksm is negligible, however because of the different 
scaling behavior of each contribution to the potential, a local minimum is possible for large 



enough h. We have plotted this case in Fig 11, again using a reference scale, L. 



IV. CONCLUSIONS 



In this work we have generalized the analysis in [TJ, demonstrating that a single brane 
embedded in a bulk spacetime can be localized and stabilized using various contributions to 
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its effective potential, specifically exploiting the non-trivial geometry and/or topology of the 
bulk manifold. We have focused on the hyperbolic horn and Euclidean cone geometries as 
they have features that are generic to more general manifolds yet still allow our calculations 
to be tractable. We have considered a brane wrapped around both manifolds, calculating 
the energy due to tension and extrinsic curvature (i.e. due to local brane geometry) and also 
due to the Casimir energy of a bulk scalar (a global effect). Both types of contributions are 
generally regional in nature, most sensitive to a finite region within the spatial manifold. The 
spaces considered here might be taken seriously for model building, as an approximation to 
capture the regional features of a more generic manifold, or perhaps simply as a toy example 
of how various forces can be used to stabilize branes embedded within non-trivial manifolds. 

We have considered both Dirichlet and Neumann boundary conditions for the bulk scalar 
at the location of the brane in both dimensions D = 2 + 1 and D — 5 + 1, i.e. without and 
with 3 Minkowski dimensions, respectively. Though we don't appear to live in D — 2 + 1, 
those systems are interesting as they illustrate how the Casimir effect changes its behavior 
with the number of dimensions, perhaps most notably how the explicit presence of poles in 
the energy depend on whether D is even or odd. In the table below we summarize our results 
on the type of minima found for each brane potential, assuming the geometric energies are 
positive (a, h > 0) and that the Casimir energy of Standard Model fields in D = 5 + 1 are 
negligible (k S m < 1). 

Case D = 2 + l D = 5 + 1 

Dirichlet Neumann Dirichlet Neumann 
Horn global none global none 
Cone global global global local 



TABLE I: Summary of brane potentials, indicating whether a global, local or no minimum is 
possible, assuming ksm is negligible. 

Interestingly, for the cone we have demonstrated that a potential minimum is possible in 
the absence of Casimir energy; the competition between purely geometrical effects (tension 
and elastic energy) are sufficient to stabilize the brane. Physically, this is the result of 
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tension energy being minimized by shrinking the volume of the brane, while the elastic 
energy is minimized by flattening it. On the cone, these two effects come in with markedly 
different dependences on the brane position, therefore an extremum is possible. For the 
horn, however, these energies have the same brane position dependence and therefore the 
same cannot be said on that geometry. 

Because the geometries we have chosen still retain a scaling symmetry, there are only 
two "special" places on each manifold, i.e. (— oo, +00) or (0, +00) on the horn and cone, 
respectively. As a result, all of the energy contributions behave monotonically with brane 
position and so are either minimized in the limit that the brane approaches one or the other. 
For this reason, we have found a necessary competition between at least two of them in order 
to achieve a stable minimum. 

In fact, this is an especially important point in the case that we live on a 3-brane. In 
those scenarios the brane may be visualized (suppressing the 3 Minkowksi dimensions) as 
a point embedded in the bulk, having only trivial geometric properties that don't result in 
position-dependent forces. With (apparently) only the Casimir effect remaining, there is 
likely no global minimum' 12 ' for the brane potential on either the horn or the cone in this 
case and such a brane could only be driven to one of the extremes of the manifold. However 
this will not generically be the case; manifolds with sufficient symmetry violation will have 
many "special" places - places where one or many of the energies are extremized, especially 
Casimir. 

Finally, we note that the discussion of regional effects of bulk manifolds is not limited 
to braneworld scenarios, alone. Perhaps the most interesting (and generic) feature of 
non-trivial regional properties of extra-dimensional manifolds is the inhomogeneity of field 
modes of the bulk. 

We would like to thank Ling-yi Huang, Harsh Mathur and Claudia de Rham for 
many useful discussions. This work is supported by grants from the U.S. Department of 
Energy and Department of Education. 



[12] However a local minimum may be possible by the competition of Casimir energies from several fields of 
varying mass and spin (see e.g. [16]). 
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Appendix A: Details of Mode Solutions 
1. Hyperbolic Horn 

a. n 7^ 



From (10) we have 

Z'Uz) ~ ZUz) + ( k 2 + 7 — n 2 e 2z ) Z n>k (z) = (Al) 



where we have defined k by u = yp 2 + A; 2 + |. With the substitutions p = \n\e z , Z(p) = 
p 1 / 2 T(p), and = — we obtain 

T pp + ^ T p-(^ + 1 ) T = (A2) 

This familiar equation is the modified Bessel equation, whose two linearly independent 
solutions are the modified Bessel functions of the first and second kind, I u (p) and K u (p), 
respectively (see, e.g. [21] )• Returning to the original variables, we have 

Z n , k (z) = e z ' 2 ( Cl I ik (Me 2 ) + c 2 K ik (|n|e 2 )) , (z > z b ) (A3) 

In this region, the Klein-Gordon norm tells us 

roc 

j dze~ z Z^zfZ^z) = 5;j (A4) 

Asymptotically, I^(x) ~ \x\^ 1 ^ 2 e x as x — > oo, regardless of p, making the solution non- 
normalizable unless we take c\ = 0. Furthermore, the modified Bessel functions I^x), K^(x) 
are monotonic for p 6 R, meaning neither K ik (|n|e 2 ) nor its derivative vanish for finite z. As 
the k -H- —k solutions are linearly related, this indicates that only real, k > are permissible. 



In fact, I u (p) and I- V {p) also form a set of linearly independent solutions to (A2), so long 
as both v ^ Z and iv Z P- 3 !. For later convenience, we use 

Z n , k (z) = e z/2 (c 1 I ik (\n\e z ) + c 2 I-i k (\n\e z ))i (z L < z < z b ) (A5) 

where the two boundary conditions at zl and z b plus the normalization condition completely 
determine C\,C2 and the allowable k. 



[13] One may worry about the possibility of k G Z, however those cases are contrived; the set of model 
parameters for which this is true form a set of measure zero. 
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For the region zl < z < z b the conditions necessary to have a zero at zl tell us (up to a 
normalization factor) 

Z n , k (z) = e z ' 2 (I^ lk (\n\e ZL )I lk (\n\e z ) - I lk (\n\e ZL )I. lk (\n\e z )) (A6) 

In the case of Dirichlet boundary condition, the zero at z b then implies the relation 

I ik (\n\e z ») = I lk (\n\e z ^) 
I- ik (\n\e*») J_ ife (|n|e^)' 

and there is a similar relation for the Neumann case. When ik G R, the monotonicity of 
the modified Bessel function implies this equation may only be satisfied in the trivial case, 
Zl — z b- All things considered, this means only real, k > are allowed. In summary, 

Z{z) = e z ' 2 [l_ ifc (Me**) l lk (\n\e z ) - I ik (|n|e 2 *) J_ ife {\n\e z )\ (z L < z < z b ) (A8) 
= e z ' 2 K lk (\n\e z ) (z b < z) 

b. n = 



From (10) we have 

Z"{z) - Z\z) + k 2 Z(z) = (A9) 
where we have defined k by uj = ^/p 2 + k 2 . The solutions are 

Z(z) = e z/2 ( Cl e ikz + c 2 e~ ikz ) (A10) 

To obviate having to delta function-normalize the modes in the region z > Zf,, we truncate 
the space at some zr > z b) taking the limit zr — > oo at the end of the calculation. This 
quantizes the modes in this region as well. After imposing zero at z^ and zr the solutions 
are 

Z ,k = e z/2 sin k(z - z L ) (z L < z < z b ) (All) 

= e z/2 sin k(z - z R ) (z b < z < z R ) (A12) 

Requiring either a Dirichlet or Neumann boundary condition in each region implies real, 
k > 0. 
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2. Euclidean Cone 

Consider the interior solutions 

Rn,k(r) = c x J v (kr) + c 2 N u (kr) (0 < r < r b ), (A13) 
for arbitrary n. Normalizability would tell us 

"drrR^rf = 1 (A14) 
If C2 7^ 0, the most dangerous part of this integral near r = is 

drrN u (kr) 2 (x J drrr~ 2u , (is ^ 0) (A15) 

oc drr (hip) 2 , (v = 0) (A16) 

Jo 

The first integral diverges because v > 1, while the second one converges. Therefore, on 
the basis of normalizibility, c 2 ^ is possible only for is = (n = 0). 

Appendix B: Calculating Eq(s) via Contour Integrals 

1. General Procedure 

We will ensure A n (k) to be a meromorphic function, having only simple zeros at the 
location of the k's in the spectrum. Near the zth zero, the derivative of its logarithm 
produces the pole (k — ki)^ 1 , plus other finite terms so that, by the Cauchy residue theorem, 
we have the relation 

£ (e + A 2 ) '^-=i-f * (* 2 + A 2 ) | In A n (k) (Bl) 

{fc} 7 

The beauty of this technique is that only the eigenfunctions need to be understood analyti- 
cally. 

Assuming no other poles/zeros in the Re > plane exist, the contour may be deformed 



(see Figure 12) by rotating counter-clockwise (k — > e m l 2 k) the upper contour, 71, while the 
lower contour, 72, is rotated clockwise (k — > e~ l7T ^ 2 k). Thus the contour goes along the 
imaginary axis from zoo to —zoo, discarding the rest of the boundary at \k\ — > 00. (This is 
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ex. k 



7i 



72 



FIG. 12: Deformed contour. 

legitimate so long as s is large enough; at the end of the calculation we analytically continue 
to s = 0). Care must be taken here as the function (k + A 2 ) 2 
a branch point at k — ±A. 
Explicitly one finds 



is multivalued and has 



1 

2ni 
1 



dk (k 2 + A 2 ) 



2ni 



ICO 

oo 



2 » lnA »« + 



m + l 



m + 1 
2 



e-™k 2 + A 2 )—~ s - (e^k 2 + A 2 ) 

cosvr(f-s) Z" 00 . o-,^-s d , . 

^ ^ / dk(k 2 -A 2 ) 2 s ^-lnA„(ifc) 

vr Ja ok 



dk (k 2 + A 2 ) 
d 



rt+l , 9 

2 * hlA "W 



9A; 



In A n (ik) 



where we have assumed that the chosen generating functions satisfy 

d d 
— \nA n {ik) = —\nA n (-ik) 



(B2) 



(B3) 



Below we shall give a brief derivation of the mode-generating functions for both the horn 
and cone geometries with the Dirichlet boundary condition. The extension to the Neumann 
condition is straightforward. 

2. Horn Generating Functions (Dirichlet, n / 0) 



Here we need to calculate 



2ni 



dk[k 2 + 



d_ 
dk 



\nA n {k) 



(B4) 
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Considering (12) we choose the generating functions 



A n {k) = -(I ik (\n\e ZL )I„ ik (\n\e z ») 



A n (k) = K ik (\n\e 



[z b < z) 



I_ ik (\n\e ZL )I ik {\n\e Zb )), (z L < z < z b ) (B5) 

(B6) 



It will be of no consequence to the actual calculation^ 14 ! but technically we had to explicitly 



include a factor of A; 1 in (B5) to avoid counting the zero mode in the spectrum (at small 



k the term in parenthesis is proportional to k). As (B5) stands it is not useful, because we 
must still perform the limit zl — > — oo. Because of the asymptotic behavior I^{x) ~ (x/2) fl 



as x — > 0, in the limit of zl — > — oo the surviving ^-dependent terms of (B5) are 



I- ik (\n\e Zb ) 
Iik(\n\e Zb ) 



on 71 (Re k > 0) 
on 72 (Re k < 0) 



(B7) 
(B8) 



When the contour is deformed, these two terms are equal and proportional to I k (\n\e Zb ) 
Thus we effectively have 



A n (ik) — I k (\n\e Zb ) 
A n (tk)=K k (\n\e z »: 



[z L <z< z b ) 
(z b < z) 



(B9) 
(BIO) 



Since the energy of each region is simply be added (as the energies of two volumes sitting 
adjacent to each other), the two generating functions combine to give 



2-ni 



dk k 2 



m + l 

IN 2 - 



d_ 

I j dk 

COS TT (y - s) 
71 



kiA n {k) 



1/2 



dk k 2 - 



(m + l) 
1\ 2 ' 



d_ 
dk 



\n[I k (\n\e Zb )K k (\n\e Zb )] (Bll) 



3. Horn Generating Functions (Dirichlet, n = 0) 

Here we will show that the sum over the /c-spectrum is ^-independent, and can therefore 
be ignored in the calculation of the Casimir energy. Before rotating the contour, we have 



[14] Since a logarithm is taken of the generating function, if we multiply it by any analytic function of k this 
only contributes an additive (though possibly infinite) constant to however this has no bearing on the 
Casimir force as it is z^-independent. 
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an integral of the form 



(B12) 



We break up the regions and, considering (All), choose the generating functions to be 
1 



Mk) = ~ k (e 



ik(zt— Zl) _ p —ik(zb—zi,y 



Mk) = ~ k (e 



ik(z R — zt,) 



— e 



-ik(z R -z b ) 



) 



[zl < z < z b ) 

(z b < Z < Zr) 



(B13) 



First we consider the zl < z < z b region, and take the limit zl — > — oo. On 71 (Im k > 0), 
Ao(k) — > e - %k { z b-z L )^-\ _ Q n ^ < A (A;) — > e^l^-^i)^- 1 . The same procedure may 

be performed for the z b < z < zr region, taking zr — > 00. Adding the two regions together, 
rotating the contour, and discarding terms that are independent of z b we are left with 



If 1 , o COS7T 

— idk k -— In A (A;) = 

27T? L afc 7T 



^ ^ - z L ) / dA; fc 



l+m— 2s 



(B14) 



The point here is that this integral is z b -independent and so not relevant for the Casimir 
effect, as claimed. 

4. Cone Generating Functions 



Considering (62) and (64), we choose the generating functions 



A n (k) = k~ v J u (kr b ) 

A n (k) = (H® (Ax^) (Ax 6 ) - (A:r out ) < 2 > (ferj)) 



(0 < r < r 6 ) (B15) 
(r 6 < r < r out ) (B16) 



The factor of k v in the top line is so that the zero mode is not technically counted 



(at small k the term in parenthesis is proportional to k u ). As it stands (B16) is not useful 
because we must still perform the limit r out — > 00. Because of the asymptotic behavior 



lim M 1 ' 2 ). 



x) ~ 




,±i(x-vTv/2-ir/£) 



TlX 



in the limit of r out — > 00 the surviving terms of (B16) are 



(kr out ) Hjp (kr b ) , on 7l (Re k > 0) 
(kr out ) (kr b ) , on 72 (Re k < 0) 



(B17) 

(B18) 
(B19) 
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When the contour is deformed, these two terms are both proportional to l v (kr b ), thus 
we effectively have 



A n (ik) = I v (kr b ) 
A n (ik) = K v (kr b ) 

Summing over the spectrum from each region we find 

(?"*) 



(0 < r < r b ) 
(r b < r < r out ) 



(B20) 
(B21) 



1 r dkk^ 9 - - ^ 



2m 



dk 



lnA n (fc) 



7T 



oo o 

dA; k l ~ 2s+m 

1/2 



0fc 



In [J„ (Ar 6 ) K v (kr b )] 
(B22) 



Appendix C: Bessel Function Asymptotics 

Using the uniform asymptotic expansion of the modified Bessel functions (see e.g. 
one may show (for the Dirichlet boundary conditions) 



(-4 + z 2 )z 2 
8(1 + z 2 f 

(-32 + 288^ 2 - 232z 4 + 13z 6 )z 2 
64(1 + z 2 f 

(-48 + 2580z 2 - 14884^ 4 + 17493z 6 - 4242,2 s + 103,2 10 )z 2 



In \I v {yz)K v {yz)\ ~ln 
where the first few fj(z) are 

U{z) 

/«(*) 



E 



fj(z) 



21), 



(Cl) 



96(1 + z 2 f 

and fj = for odd j. Likewise, (for the Neumann conditions) one may show 



(C2) 
(C3) 
(C4) 



(4 - 3z 2 )z 2 
8(1 + z 2 f 

(32 - 32(k 2 + 328z 4 - 27z & )z 2 
64(1 + z 2 f 

(48 - 2652^ 2 + 1618(k 4 - 207992 6 + 5652^ 8 - 162^ 10 ),2 2 



hi[K' v {uz)l' v {vz)\^ 

where the first few gj(z) are 

92{z) = 

9a{z) = 

%{z) = 
and gj = for odd j. 



i Vz^Ti 

~2v~ J 2 



E 



9j(z) 
pi 



(C5) 



96(1 + z 



2\9 



(C6) 
(C7) 
(G8) 
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Appendix D: Details of 2 + 1 Horn Calculation 



1. Ef-(s) 



Renormalization of Brane Tension: We see in equation (35) that a pole remains' 15 ', 
specifically one that is z^-dependent: 

9 1 

— ^_ e -*6 x _ ( D1 ) 

256 s y 1 

The Casimir force cannot be infinite, so clearly we must identify a physical parameter to 
absorb this divergence. Notice that this term is proportional to the volume of the brane 

i-2-k 

KraneOC / d0y/fr\ = = 2716^ (D2) 

Jo 

This suggests that the brane tension is a renormalizable quantity that will suffice. Let us 
consider the energy due to a bare brane tension, a: 

r-2jT 



E ten = / d9^yi\a = 2ne- zb x a 
Jo 



(D3) 



For convenience we define the renormalized brane tension using not only terms from Eq 1v , 
but also the finite contributions to E which are proportional to the brane volume (in the 
limit z b — > — oo): 



E u , x = 2k< = 2k< : 7 ^ ^ (~ + 7 + y s + fini1r ) ) [ D U 



■)) 



With this definition 



lim E{z h ) ~ E tcn = 2ne zb x a ren , (D5) 



which is now finite and independent of \x. 



2. ££ em 



From (36) we have 



n=l 



1 + 16x 2 b + 49x£ + 6Axf, . (\ 
4(l+4^) 5/2 



+ ^6 + 



128' 



(D6) 



[15] See also Sec. E 2 a for an independent check using the heat kernel technique. 
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and remind the reader that x b = \n\e Zb . 



Limit z b — > oo: Since x b ^> 1 for all n, we can expand in x b 1 so that 



OO 



lim E™ m ~-Y —xf = -— C(3)e- 3 ^ (D7) 
zb^oo ^ 1024 b 1024 s v ; 

n=l 



Limit —7- — oo: We can approximate the sum as an integral by using the Euler-Maclaurin 
formula which, to lowest order, says (see e.g. |25j) 

f(An) « - / dxf(x) + (f(A) + f (oo)) (D8) 

n=l A jA Z 



We identify A e* 6 so that f(x b ) is the summand of (D6), thus 



lim E T em ~ e - z » [ dxf(x) + lf(e Zb ) 

Z{,->-oo J eZ6 2 



7+ 181n[4e 2 H z 1 /ri/ 2z N 
~ e +o+ 0{e 2zb ) 



256 



(D9) 



3. 



From (38) we have 



(D10) 



(* 2 -i)' 
where 

£> fc (a; 6 ) = (in [K k (x b ) I k {x b )\ - In [K k (x b ) I k (a*)]) (Dll) 

Limit ^ — 7- oo: Again, it will always be true that x b ^> 1. Under this condition, the uniform 
asymptotic expansions for the modified Bessel functions are a good approximation because 
beyond the leading logarithmic term they are expansions in inverse powers of x b . Therefore 
by going to the next highest order in x^ 1 we find (see Appendix O 



Uk [Xb) ~ 64(sg + F)^ 
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This can be directly integrated and, after performing another expansion in x b 1 



b 

OG 



E*« « 



37tt r 
16384 Xb 



n=l 

S77T 

-^C(3) 6 - (D13) 



Limit ^ —7- — oo: Once again, we will use the lowest order Euler-Maclaurin expansion. In 
this limit one may show that 

limDJa;) ~ 0(x 2 ,x 2k ) 

x—>0 

(D14) 

depending on whether k is less than or greater than 1. This implies that to good approxi- 
mation 

e~ Zb f°° „ k 
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E o n -—I dk — " / dx D k (x). (D15) 



h/2 (k*-lY 

In the above expression we have lowered the lower limit of integration to zero, the error 
being roughly 



re o ne o 

e~ Zb dx D k (x) ~ e~ Zb / dx 0(x 2 , x 2k ) ~ 0(e 2zb , e 2kZb ) (D16) 
Jo Jo 

which is negligible as Zb — > — oo. We have thus learned than in this limit E^ n oc e~ Zb , so it 
remains to calculate this proportionality factor. Numerically computing the double integral 



in (D15), we find 

E Rn « 2.4 x 1(T 3 x e~ Zb (D17) 

Appendix E: Heat Kernel Analysis 

In this appendix we give a brief review of the heat kernel technique and apply it to 
confirm the poles encountered in D = 2 + 1, providing a partial check on our calculations. 

1. Heat Kernel Analysis of a General Eq(s) 

Following the nice review by Vassilevich [26] and book by Fursaev and Vassilevich [T9] . 
if we use a Mellin transform to write 

POO y.s-3/2 

u[~ 2s = / dt-) ^e-^ 2 (El) 

Jo 1 [s - 2J 
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then the zeta-function regularized Casimir energy (|4]) can be recast as 



2s roc y.s-3/2 

= ^r/ *t— r-K® (E2) 

2 JO r ( S "2) 

where K(t) is the heat kernel trace, given 

by[16] 

K(t) = e ~ twf = tE ^ a p ( E3 ) 

i p=0 

where n is the number of spatial dimensions and the heat kernel coefficients, de- 
termined by integrals over local geometric scalars and also depend on the field type and 



boundary conditions. Considering (E3), in the limit s — > the divergent terms of (E2) are 
those in which p < n + 1. Analytically continuing s — > gives a finite answer for all terms 
except when p — n + 1; this is an inescapable pole. The divergence comes from the t — > 



part of the integral, so we break the integral (E2) up into two regions, [0, 1] and [1, oo). The 



later converges, while the former contains the pole: 

fx 2s f 1 t 3 - 1 fx 2s t 



o r( s -§) n+1 4^ 



s 



1 







+ O(s ) 



=On+l- ( E4 ) 



A^/n s 

2. Heat Kernel Technique Applied to the D = 2 + 1 Geometries 

To verify the poles in the Casimir energy encountered for both the horn and cone geome- 
tries in D = 2 + 1 (n = 2) we require the coefficient 03, for which we borrow the results 
from , suited for our problem and conventions. To apply them, we identify the variables 
from that work, (/, E, S, U + , Il_) as / = 1, E — 0, and 5 = 0, with n + and Il_ depending 
on the choice of boundary condition. With these results in hand, we have 

a 3 = -J-= [ d^xy^l (16 X R ~ SRanan + (1311+ - 7I1_) K 2 + (2n+ + 10n_) K ab K ab ) 
768V7T J dM 

(E5) 



[16] Technically the summand should be e tXi , where the Xi are the eigenvalues of the generalized Laplace- 
Beltrami operator. Since we use a simple massless Klein-Gordon field, ujf = A.;. 
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where 7 M „ is the induced metric on the boundary, x = n+ — the Riemann tensor is 
defined^ by [V M , VJ A a = R a Xllu A x , and K ab is the extrinsic curvature tensor of the 
brane. Latin (or hatted) indices designate orthonormal coordinates, and repeated indices 
imply a summation. As we have split up both of the geometries into two regions, we take care 
to calculate 03 one each side of the boundary separately, however in the cases we consider 
the answers are identical. 



a. Hyperbolic Horn 

Here the determinant of the induced metric on the brane is given by I7I = e~ 2zb , where 
we have set the horn curvature scale, = 1. The relevant geometric objects are 

R = -2 (E6) 

Rom = -1 (E7) 
K ab = ±S aS 5 b§ (E8) 

where the sign of K a i is positive (negative) if the normal to the boundary, points in the 
positive (negative) z direction. 

Dirichlet Condition 

Here we identify n + = 0, n_ = 1. For either region we find 



Considering (E4), we confirm that the entire pole part of the energy is 



E » M = < E10 > 



[17] There is a difference of a minus sign between the Riemann tensor as defined in [26j and here. 
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Neumann Condition 



Here II + = 1, II_ = 0. For either region we find 

a * = -^^ e ~ Zb ( EU ) 

implying 

EZ(.s) = A e ^i ( E12 ) 



6. Euclidean Cone 

Here, the determinant of the induced metric on the brane is given by |7| = r\. As the 
geometry is Euclidean, only the extrinsic curvature tensor is non-zero: 

K ab = ±-5 a9 5 be (E13) 
n 

where the sign of K ab is positive (negative) if the normal to the boundary, n M points in the 
positive (negative) r direction. 



Dirichlet Condition 



Here n + = 0, n_ = 1. For either region we find 

= (l-*)y/i 

128r fc V ; 



E%(s) = -^J 1 - (E15) 



Considering (|E4|), we confirm the entire pole part of the energy is 

(1-5)1 
256r b s 

Neumann Condition 



Here n + = 1, n = 0. For either region we find 

y-5) 

128r b 



a 3 = 5 ' 1 -^>^ (E16) 



implying 



< W = - 5 li^)i (E17) 
V ; 256r 6 s V ; 
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